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ABSTRACT 

Hubble Space Telescope images of SCR radio galaxies at redshifts 0.6 < z < 1.8 
have shown a remarkable variety of structures, generally aligned along the radio axis, 
indicating that the radio source strongly influences the optical appearance of these 
galaxies. In this paper we investigate the host galaxies underlying this aligned emission, 
combining the HST data with ground-based infrared images. 

An investigation of the spectral energy distributions of the galaxies shows that the 
contribution of the aligned blue component to the K-band light is generally small 
(~ 10%). The radial intensity profiles of the galaxies are well matched at radii < 35 kpc 
by de Vaucouleurs' law, demonstrating that the K-band light is dominated by that of 
an elliptical galaxy. There is no evidence for a nuclear point source, in addition to the 
de Vaucouleurs profile, with a contribution > 15% of the total K-band flux density, 
except in two cases, 3C22 and 3C41. We conclude that the K-band emission of the 
distant 3CR galaxies is dominated by starlight. The magnitudes, colours and location 
of the distant 3CR galaxies on the projected fundamental plane indicate that their 
stellar populations formed at high redshift and have since been evolving passively. 

Large characteristic radii are derived for the 3CR galaxies, indicating that they must 
be highly evolved dynamically, even at a redshift of one. At radii larger than ^ 35 kpc, 
a combined galaxy profile clearly shows an excess of emission as compared with de 
Vaucouleurs' law, indicating that at least some of the galaxies possess cD-type halos. 
This supports other independent evidence for the hypothesis that the distant 3CR 
galaxies lie in moderately rich (proto-)clusters. Since the nearby FRII galaxies in 
the 3CR catalogue lie in more diffuse environments and do not possess cD halos, the 
galactic environments of the 3CR galaxies must change with redshift. The K— z relation 
of the 3CR galaxies cannot, therefore, be interpreted using a standard 'closed-box, 
passively evolving stellar population' model, whereby the galaxies that host distant 
3CR sources will evolve into the galaxies that host nearby 3CR FRII sources. 

At redshifts z ~ 1, the absolute K-magnitudes of the stellar populations of the 3CR 
galaxies are brighter than those of the lower radio power 6C galaxies, indicating that 
the 3CR galaxies contain a greater mass of stars; this is consistent with them lying 
towards the centres of clusters. Powerful high redshift radio galaxies possess radio 
beams whose kinetic power is close to the Eddington limiting luminosity of a central 
supermassive black hole. Since the mass of the black hole is likely to scale in proportion 
to the mass of the host galaxy, the 3CR galaxies will contain more massive central 
engines than the 6C galaxies, which accounts for their more powerful radio emission. 
At redshifts z < 0.6, the beam power of the radio sources is limited by the availability 
of fuel for the central engine rather than the black hole mass, and so no correlation is 
expected between the radio power and the mass of the host galaxy. 
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1 INTRODUCTION 

The 3 CR sample of radio sources defined by Laing et al. 
(1983) consists of the brightest extragalactic radio sources 
in the northern sky, selected at 178 MHz, and contains radio 
galaxies and radio quasars with redshifts up to 2 ~ 2. The 
host galaxies of the low redshift sources in the sample are 
identified with giant elliptical galaxies containing old stellar 
populations and, if the host galaxies of the high redshift 



sources rSTB' 

can be used Lu sLud^ the evuluLiuu of Lhis class uf galaxy 
with cosmic epoch and consequently to constrain models of 
galaxy formation and evolution. 



infrared emission may contain light directly from the central 
AGN; (iii) the observation that infrared emission lines such 
as [SIII] 9532 can contribute a sign ificant percentage of the 
K-band flux (Rawlings et al. 1991). 

In addition, the K— z relation of the 6C radio galaxies, 
which are roughly a factor of five lower in radio luminosity 
than the 3CR sources, tends to track that of the SCR galax- 
ies at low redshift, but at higher redshift the 6C gala xies tend 
to be fainter, lying closer to the ' no evolution' line (Bales & 



also giam elhpLicals, .hen- sLellar populadons Rawhngs 199| [Bales et al. 1997|). This result is consistent 



Lilly and Longair (1982; 1984) obtained infrared K- 
magnitudes for an almost complete sample of 83 3CR galax- 
ies with redshifts < z < 1.6, and constructed the K- 
magnitude vs redshift relation for these objects. The re- 
sulting relation showed remarkably little scatter and was 
interpreted as indicating that the 3CR host galaxies at red- 
shift z ^ 1 are indeed giant elliptical galaxies. Lilly & Lon- 



gair (1984) also showed that, unless the deceleration pa- 
rameter were as large as qo ~ 3.5, the shape of the K—z 
relation would not be consistent with non-evolving stellar 
populations, but that at least passive evolution is required. 
The K—z relation suggests that the host galaxies must have 
formed at large redshift, ^ 3; this age is consistent with 
the red colours of the infrared emission of some of these 



galaxies ( Lilly 1989 ). Support for such an early formation 
epoch comes from a very deep spectrum of 3C65, at redshift 
z = 1.176, which is well -matched using a ste llar population 
of at least 4 Gyr in age ( Stockton et al. 1995 ), and also from 
the recent results of Dunlop et al. (1996). 

In 1987, McCarthy et al. and Chambers et al. discov- 
ered that the optical emission of these powerful radio galax- 
ies tends to be aligned along the axis of the radio source. 
Many models have been proposed to account for this align- 
ment effect, the most promising being massive star forma- 
tion induced by the passage of the radio jets (e.g. Rees 1989) , 
scattering of light from an obscured active galactic nucleus 
(AGN) by electrons or by dust (Cimatti et al. 1996,1997, 
Dey et al. 1996, and references therein) and nebu lar con- 
tinuum emi ssion from warm line-emitting regions ( Dickson 
et al. 1995). This complicates the use of these galaxies as 



cosmological probes. 

A number of observations have also suggested that a 
proportion of the K-band light of these galaxies may not 
be associated with starlight: (i) the discovery of a weaker 
alignment effect at near-infrared wavelengths (e.g. Eisen- 
hardt and Chokshi 1990; Rigler et al. 1992; Dunlop & Pea- 
cock 1993), which may be produced by the long wavelength 
tail of a fiat spectrum component responsible for the opti- 
cal alignment effect; (ii) the detection of broad Ha emission 
from the d i stant radio galaxy 3C2 2 at z = 0.938 (Rawlings 



with the suggestion of Yates et al. (1986) that there is an in- 
trinsic correlation between the absolute infrared magnitude 
and the radio luminosi ty of powerful radio galaxies at large 
redshifts. Eales et al. (1997) have argued that a proportion 



et al. 1995; Economou et al. 1995) and a reported detection 
of an unresolved centr al component in the K-band emission 
of 3C65 at 2 = 1.176 ( Lacy et al. 1995|) , suggesting that the 



of the K-band emission of the 3CR galaxies may be directly 
or indirectly associated with AGN activity, the K—z relation 
being caused, in part, by the correlation of radio luminosity 
with redshift in the fiux-limited 3CR sample. 

In order to study the astrophysics of powerful distant 
radio galaxies, we have selected an almost complete sample 
of 28 sources from the 3CR catalogue, in the redshift range 
0.6 < z < 1.8, for observation by the Hubble Space Tele- 
scope (HST). In a previous paper (Best et al. 1997; hereafter 
Paper I) we presented the results of observations of these 
sources at optical wavelengths using the HST, at 8.4 GHz 
using the Very Large Array radio interferometer (VLA), and 
at 1.2 and 2.2 ^m using the IRCAM3 array on the United 
Kingdom InfraRed Telescope (UKIRT). The HST observa- 
tions showed that the optical morphologies of the major- 
ity of the radio galaxies bear little resemblance to giant el- 
liptical galaxies, being dominated instead by high surface 
brightness structures elongated along the axes of the dou- 
ble radio sources (Paper I) . In interpreting the HST images, 
however, it should be appreciated that standard elliptical 
galaxies possess very red colours and, at a redshift of one, 
have low surface brightnesses at optical (rest-frame ultravi- 
olet) wavelengths (e.g. Giavalisco et al. 1996). Only regions 
of enhanced ultraviolet emission appear prominently in the 
images. 

We defer an analysis of the variety of optical structures 
and of the alignment effect to the third paper in this series 
(Best et al. in prep.; Paper III). Here, we restrict the in- 
vestigation to the old stellar populations of the host radio 
galaxies. The outline of this paper is as follows. In Section ^ 
we describe the data reduction used in the current analysis 
which has not already been discussed in Paper I. In Sec- 
tion ^ we address the problem of estimating the fraction of 
the K-band light which is associated with a flat spectrum 
aligned component. In Section ^ we investigate the radial in- 
tensity proflles of the galaxies. The relationship between the 
half-light radii and the surface brightnesses of the galaxies 
is compared with the locus of low redshift elliptical galaxies 
and brightest cluster galaxies in Section ^. In Section ^ we 
derive an improved K—z relation for the narrow line radio 
galaxies. The implications of our results for the formation 
and evolution of massive galaxies and for the cosmic evolu- 
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Table 1. Measured properties of the radio galaxies (see Paper 
I). Column 1 contains the 3C catalogue name of the source, with 
its redshift in column 2. The K-magnitude measured from our 
UKIRT images, is given in column 3, with the error in column 
4. This was generally measured through a 9 arcsec aperture (for 
those galaxies indicated by an asterisk, a 5 arcsec diameter aper- 
ture was used to exclude nearby companions; see Paper I), and 
has been corrected for galactic extinction using the extinction 
maps of Burstein and Heiles (1982). For 3C368 (marked **) the 
K magnitude quoted is after subtraction of the M— star (see Paper 
I). The ellipticity of the K-band image is given in column 5 and 
the position angle offset from the radio axis in column 6. In col- 
umn 7 we tabulate the Kcorr magnitudes obtained after scaling to 
a 9 arcsec diameter aperture and correcting for any fiat spectrum 
and point source contributions (see text). 



tion of the radio source population are discussed in Section [7[ 
and we summarise our conclusions in Section ^. 

Throughout the paper, t he GI SSEL stellar synthesis 
codes of Br uzual and Chariot (1995) are used. We adopt a 
Scalo ( |l986| ) initial mass function with upper and lower mass 
cut-offs of 65Mq and O.IMq respectively, and solar metallic- 
ity. Except where otherwise stated, a deceleration parameter 
go = 0.5 and a Hubble constant Ho = 50kms~^ Mpc~^ are 
adopted. 



2 THE DATA 



A variety of infrared properties of the galaxies were mea- 
sured directly from the infrared images and are presented 
in Table |l| (see also the data tables presented in Paper I). 
Table ^ contains properties of the galaxies which are derived 
in the text. 



The photometric magnitudes of the galaxies were mea- 
sured at a variety of wavelengths (see Paper I), and were con- 
verted into broad-band flux densities for the investigation of 
the spectral energy distributions (SEDs) of the galaxies in 
Section ^. Most of the optical flux densities are influenced to 
some extent by emission lines which lie within the passband 
of the filter. In Paper I, we calculated the percentage of the 
flux density that is associated with line emission for each 
of the HST observations; these values are used to produce 
'true' broad-band continuum flux densities. In addition, the 
contributions of the [OIII] 4959,5007 and H/3 emission lines 
to the J-band flux densities of those galaxies at redshifts 
z > 1.1 have been estimated from the [Oil] 3727 line flux, 
using line ratios taken from the c ombined radio galaxy spec- 
trum derived by McCarthy ( |l993|) , and have been subtracted 
from the J-band magnitudes. No correction has been made 
for line contamination in the K waveband, although this 
should produce little error since strong lines such as Hq do 
not enter the K-waveband until a redshift of 2 ~ 1.9, which 
is greater than that of any of the sources in the sample. The 
only emission line of any signiflcance which may be present 
in the K-band is the [SIII] 9532 line. Rawlings et al. ( |l99l| ) 
detected a signiflcant (~ 15%) contribution to the K-band 
flux density of 3C356 in this line, but in none of the other 
galaxies they observed was its contribution greater than ~ 5 
to 10%, and often the line was not detected. In addition, 
this line only falls in the K-band for the eight galaxies with 
redshifts 1.05 < z < 1.55, and these galaxies do not have 
signiflcantly brighter absolute K-magnitudes. 



Another factor which influences the measured flux den- 
sities is reddening by dust intrinsic to the galaxies. Since no 
accurate estimates have yet been made of the mass, tem- 
perature and distribution of dust throughout these galax- 
ies, taking account of dust extinction is problematic. It has 
a negligible effect at the long wavelengths of the K-band 
emission, with which most of this pa per i s concerned, but it 
will affect the SED fitting in Section 



3.1 



The distribution of dust within elliptical galaxies can 
be divided into two components. First, there is dust spread 
throughout the galaxies: this affects the colours of the galax- 
ies but has little effect upon their morphologies. This dust is 
therefore impossible to detect in our observations, and can- 
not be properly compensated for. We discuss this extended 
dust briefly when modelling the SED in Section 3.1. The 
second dust component is that whic h lies in disks near the 
centre of the galaxies; de Koff et al. (1996) showed that be- 
tween 30 and 40% of low redshift radio galaxies possess clear 
dust lanes. The most prominent example of this phenomenon 
in the galaxies in our sample is the large dust lane running 
across the central regions of 3C324. Dickinson et al. (199f ) 
showed that an extinction E(B— V)~ 0.3 is required in the 
central regions of this galaxy to obscure the galaxy nucleus; 
at the wavelengths of the HST observations, this extinction 
coefficient corresponds to removing about 70% of the flux 
density from the regions underlying the dust lane (e.g. de 
Vaucouleurs & Buta, 1983). For 3C324, this correction re- 
sults in a 19% increase in the total flux density through the 
F702W filter, and 16% through the F791W filter. 

Four other galaxies, 3C68.2, 3C252, 3C266 and 3C356, 
also show reasonably convincing evidence for dust lanes, al- 
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though not to the extent of 3C324Q|. The much more com- 
pact dust lanes in these sources, relative to that of 3C324, 
make accurate measurements of the extinction coefficient 
impractical, and so we adopt a value of half of that of 3C324. 
Repeating the above procedure for these galaxies provides 
best estimates of the optical flux density corrections of, in 
order of increasing wavelength for each source: 3C68.2 — 
7%; 3C252 — 7%, 4%; 3C266 — 10%, 5%, 4%; 3C356 — 
6%, 4%. These corrections have been incorporated into the 
flux density measurements used in the next section. 



3 THE INFRARED ALIGNED COMPONENT 

3.1 Fits to the spectral energy distributions 

The infrared observations presented in Paper I confirm the 
result that the K-band emission of the powerful radio galax- 
ies in our sample has a tendency to be aligned along the radio 



1992 



Dun- 



axis ( Eisenhardt fc Chokshi 199C ; Rigler et al 
lop fc Peacock 199(^ ) but it is a significantly weaker effect 



than iui thc optical waveband. Thia aligned infrared cmioaion 
may bd the long wavelength tail of the flat spectrum com — 
ponent responsible for the alignment effect in the optical 
waveband, or may indicate that the old stellar population 
itself is aligned with the radio source, as in the models of 



Bales ( |1992| ) and West ( |l99l| ; |99|). 

A relatively simple procedure has been adopted for es- 
timating the contribution of the optically aligned compo- 
nent to the infrared alignment. For each galaxy an SED is 
produced which best matches the observed broad-band flux 
densities, using the sum of two components. The flrst compo- 
nent is an old stellar population, whose SED is derived using 
t he G ISSEL stellar synthesis codes of Bruzual and Chariot 
(1993). We assume that the stars formed in a 1 Gyr burst 
which occurred at redshift z — 10, and the stellar population 
then evolved passively until the observed epoch; negligible 
changes result for the SED at redshift z ~ 1 if a formation 
redshift of 5 or 50, or an exponentially decreasing star forma- 
tion rate are used instead. The second component represents 
the aligned emission; the spectral shape of this component 
depends upon the physical cause of the alignment effect. 

For electron scattering, a power-law quasar spectrum is 
expected, oc , where the value of a is generally about 
—0.7 to —1 for quasars (e.g. Baker & Hunstead 1996). Due to 
the much lower scattering efficiency of electrons compared to 
dust (e.g. di Serego Alighieri et al. 1994), if electron scatter- 
ing is to dominate, then the mass of extended dust must 
be small, and extinction by the extended dust will have 
Uttle effect on the spectral shape. If scattering by dust is 
entirely responsible for the excess ultraviolet fiux then, as- 
suming the dust to be evenly distributed throughout the 
galaxies, a few times lO^M© of dust are required, (e.g. Di 



* Note that there is marginal evidence for weaker dust lanes in 
the central regions of other galaxies in the sample; however, these 
are too weak for reliable estimates of their effect to be made and, 
in any case, the increase in flux density would be < 5%. De Koff 
et al. (L996) suggested that the frequency of dust lanes in the 



SCR radio galaxies with redshifts 0.1<2<0.5is invariant with 
redshift, corresponding to about 30 to 40% of all the sources. Our 
images would be consistent with a similar percentage of dust lanes 
being present in the high redshift SCR galaxies. 



Serego Alighieri 1994). Light scattered by this dust will have 
a spectral s hape somewhat blue r than the incident q uasar 
power law (Calzetti et al. 1995; [Cimatti et al. 1996 ), but 
dust outside the scattering regions will redden the scattered 
emission, giving it a similar spectral shape to electron scat- 
tering. A young starburst has a spectral shape which is rel- 
atively fiat in whilst nebular continuum emission has a 
much redder spectral shape below about rest-frame 3600A 
but provides little K-band intensity. Each of these effects is 
likely to be present at some level, and so the overall SED of 
the aligned component will be some combination of these. 
For illustrative purposes. We will assume that one third of 
the aligned emission at rest-frame 3600A is associated with 
each of the three processes of scattering, star-formation, and 
nebular continuum; this is well within the range of previous 
estimates (e.g. Dickson et al. 1995, Cimatti et al. 1997, Tad- 
hunter et al. 1996). This combination can be approximated 
from optical wavelengths through to K-band by a single 
spectral shape, cx r/"" *, which is in close agreement with 
the fi, oc u"^'^ rest-frame ultraviolet spectral s hape deter 



mine d for the powerful radio galaxy 4C41.17 (Dey et al 
1997). This, fv oc v'" *, is the spectral shape we adopt for 



the aligned emission, and hereafter we refer to this compo- 
nent as the 'fiat spectrum component'. 

For each galaxy, the old stellar population and flat spec- 
trum components were weighted to produce the combined 
spectrum which best matched the broad-band fiux densities 
measured in the UKIRT and HST observations. These fits 
are shown in Figure ^. Such fits were not made for 3C437 
and 3C470, for which only two broad band fiux densities 
were available, nor for 3C22 and 3C41 which both possess 
significant nuclear components in the infrared images (see 
Section 1^. 

The fits shown in Figure |l| are, in the majority of 
cases, rather good, indicating that the simple two compo- 
nent model of the SED can provide a good representation 
of the spectra of these galaxies. The mass of stars within 
the galaxy and the fiux density of the fiat spectrum compo- 
nent are listed in Table S. Notice that the stellar masses are 
typically a few times 10 ^ Mq which, given that the mass- 
to-light ratios of such galaxies are usually well in excess of 
1, demonstrates that these galaxies are significantly more 
massive than M* . The contributions to the K-band fiux of 
the two components of the fit can be compared to provide an 
estimate of the fraction of the emission in this band which is 
associated with the aligned component. The values derived 
from our best fitting models are tabulated in Table ^ and 
range from ~ 1% in the passive sources 3C65 and 3C337, 
up to ~ 31% in the case of 3C368, with an mean value of 
about 10%. 



3.2 The origin of the aligned infrared emission 

To quantify the degree of alignment, we define the alignment 
strength as: 



1 



45 y 



where A6 is the difference in position angle between the K- 
band emission and the radio axis, measured in degrees, and 
e is the ellipticity of the K-band image, defined as: 
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Figure 1. SED fits to tiie broad band flux densities of the SCR galaxies using an old stellar population and an aligned component (see 
text and Table H for details). 



where p and q are dimensions along the major and minor 
axes respectively of the infrared image of the radio galaxy. 
The alignment strength therefore ranges from +1 to — 1: 
those galaxies in which the major axis of the infrared emis- 
sion is oriented within 45° of the radio axis have a positive 
alignment strength, whilst those which are misaligned have a 
negative value; a value of zero corresponds either to a galaxy 
whose major axis is misaligned by 45° from the radio axis 
or to a galaxy which is circularly symmetric. This is a bet- 
ter statistic than considering Ad alone since it results in a 
higher value for a source which is highly elongated but mis- 
aligned by, say, 15° from the axis of the radio emission than 



for a source which is virtually symmetric but whose major 
axis lies directly along the radio axis. The values of the K- 
band alignment strengths are tabulated in Table ^ and are 
plotted as a histogram in Figure ^. There is a clear tendency 
for the galaxies to show some alignment in the K waveband, 
with nearly 70% possessing positive alignment strengths, a 
result significant at greater than the 95% confidence level. 
The mean value of the alignment strength is 0.05. 

It is interesting to compare the K-band alignment 
strength with the percentage flat spectrum contribution 
in the K-band (Figure Of the five galaxies with the 
largest K-band flat spectrum components, four are amongst 
the seven sources with high K-band alignment strengths, 
as > 0.10, including all three in the highest bin in Figure H 
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Figure 1. cont. 



(see also Table g). A Spearmann rank test indicates that 
these parameters are correlated, at a significance level of 
about 97.5%. If the four sources, 3C68.2, 3C252, 3C266, and 
3C368, with particularly high values of both as and the flat 
spectrum contribution are excluded, then for the remaining 
24 galaxies the average alignment strength is (ss) = 0.02, 
different from the hypothesis of no alignment at only low 
significance. 

The correlation displayed in Figure |^ indicates that the 
long wavelength tail of the active blue emission responsi- 
ble for the alignment of the optical images is an important 
source of the K-band alignment effect. Scatter in this rela- 
tion will arise from the errors introduced in determining the 
two parameters from the data and, more importantly, from 
any intrinsic scatter in the alignment strength, such as that 



suggested by the negative tail in Figure ^ An estimate of the 
underlying scatter can be gauged from the range of align- 
ment strengths in Figure ^ of those galaxies with only small 
flat spectrum components, suggesting a standard deviation 
of between 0.05 and 0.1. If the parent galaxies are giant el- 
liptical galaxies, this can naturally be explained in terms of 
the radio axis being randomly orientated with respect to the 
major axis of the galaxy. 



4 RADIAL LIGHT PROFILES 

Our observations can be used to test whether the K-band 
light distributions of the galaxies have the profiles of el- 
liptical galaxies, by attempting to fit de Vaucouleurs' law. 
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Table 2. Properties of the radio galaxies derived in the various sections of the paper. Co lumn 1 contains the 3C catalogue name of the 
source. Column 2 lists the stellar mass of the galaxy derived from the fits in Section and column 3 the flux density at rest— frame 
3000A of the 'flat spectrum' component derived from the same flt. The values for 3C437 and 3C470 (marked with an asterisk) are only 
estimates since only two broad band flux densities were available. This analysis was not possible for 3C22 and 3C41 because of a nuclear 
contribution to their flux densities. For the remaining 24 galaxies, the approximate percentage of the K-band light that is associ ated 
with a flat spectrum aligned component is given in column 4. The K-band alignment strength of each galaxy, as deflned in Section ^.2[ 
is listed in column 5. Column 6 gives the best-fltting de Vaucouleurs radius for each galaxy (see Section W, where the errors involved are 
also discussed). Such flts were not possible for 3C247 and 3C368, nor for the flve highest redshift sources (see text). The percentage of 
K-band light associated with an unresolved nuclear emission source in the best-fltting model is given in Column 7. The value after a 
correction factor has been applied (see Section ^ for further details) is given in Column 8, and the 90% confldence limits for this value 
are given in column 9. The measured mean surface brightness within the de Vaucouleurs radius in the K-band is given in column 10, 
and the modelled V-band surface brightness for a passive evolution model is listed in column 11 (see Section ^ for more details). This 
analysis was not carried out for 3C22 and 3C41 because of the presence of a nuclear contribution to the K-band flux. 

7(r) oc exp [—7.67 (r/rc)^^^] , to their radial surface bright- 
ness profiles. The HST images allow such plots to be made 
without the need to account for the effects of seeing, and 
would allow an excellent determination of the radial profiles 
to be made, were it not for the fact that at optical wave- 
lengths the emission from the majority of the galaxies is 
dominated by the aligned blue component and so this anal- 
ysis cannot be carried out. Eight of the galaxies, however, 
do not show strong ultraviolet emission. The HST images 
show that they possess almost symmetrical morphologies 
and only weak aligned components are required to fit their 
SEDs (< 10% of the total flux density in the HST image). 
The radial intensity profiles of the optical emission of these 
galaxies were measured, nearby companions objects being 
removed and replaced with the average background intensity 
at that distance from the centre of the galaxy. No compan- 
ion subtended an angle greater than 20 degrees about the 
centre of the galaxy, and so any errors in this replacement 




-0.4 -0.2 0.2 

Alignment Strength 



0.4 



Figure 2. A histogram of the K-band alignment strengths. Val- 
ues lying on a boundary are included in the lower bin. 
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Figure 4. De Vaucouleurs fits to the radial intensity profiles of the HST images of eight SCR radio galaxies which do not show a 
significant active ultraviolet component. The characteristic radius of each, determined from the gradient of the best fitting straight line, 
is given. 
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Figure 3. A plot of the K-band alignment strength against 
the percentage contribution at K of the fiat spectrum component 
required to match the SED of the galaxy. 



will be small. The radial profiles are displayed in Figure ^ 
They show that a de Vaucouleurs profile provides an excel- 
lent fit to the observed radial surface brightness in six of the 



eight cases, providing some of the strongest evidence to date 
that these high redshift galaxies are indeed giant elliptical 
galaxies. The derived values of rc are listed in Table ^ and 
have a typical accuracy of about 15%. The cases of 3C22 
and 3C41, the two sources for which a good de Vaucouleurs 
fit is not obtained at small radii, are discussed below. 

The characteristic radii, rc, obtained from these fits 
were then used to investigate the infrared K-band radial in- 
tensity profiles. Where companion galaxies subtended less 
than about 45 degrees with respect to the centre of the 
galaxy, these were removed as in the HST images, and re- 
placed by the average of the background pixels at that point. 
In two cases, 3C34 and 3C441, nearby companion galax- 
ies subtended an angle > 45° about the centre of the 3CR 
galaxy, and this substitution could lead to large errors in 
the radial profiles. In these cases, the radial intensity profile 
out to a radius of about one arcsec was determined using the 
normal method, but at larger radii only the 2 or 3 quadrants 
around the source which did not contain the companion were 
used, and the emission from these was then scaled to derive 
the average intensities. The two methods matched well at 
the boundary at 1 arcsec radius. 
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Figure 5. Fits to tlie radial intensity profiles of tlie K— band images of the eight 3CR galaxies shown in Figure using the sum of an 
unresolved point source (dash-dot line) and a de Vaucouleurs profile with the characteristic radius determined from the HST images 
(dashed line). For each of the profiles, the effect of seeing has been taken into account. The sum of the two components is indicated by 
the solid line. In many cases the best fit does not involve a point source component, and so only the solid line is shown. 



The infrared radial intensity profiles of these eight 
galaxies were modelled using a two-component fit, involv- 
ing a de Vaucouleurs profile with the characteristic radius 
determined from the fits to the HST images, and a point 
source of emission; each profile was convolved with a Gaus- 
sian profile matching the seeing of the infrared observations. 
The seeing, which was typically between 1 and 1.2 arcsec, 
was measured using objects on the infrared images of the ra- 
dio galaxy which appeared unresolved in the HST data. The 
fitting procedure involved minimising the sum of the squares 
of the offsets of the combined model profile from the data 
points within the inner 3 arcsecond radius, weighting each 
data point inversely by its error. Tests showed that the pro- 
cedure converged to the same minima from a complete range 
of initial values. 

The resulting fits are shown in Figure ^: the dashed line 
shows the de Vaucouleurs profile, the dotted line shows the 
radial profile of the point source, and the solid line shows 
the total intensity profile. In many cases, the point source 
component is zero, and so only the solid line is plotted. In 
each case a good match is obtained. The percentage of the 



total K-band emission of each source associated with a point 
source for the best fitting combination is given in Table ^. 

This analysis indicates the reason for the poorness of fit 
of the de Vaucouleurs profiles to the optical radial intensity 
profiles of 3C22 and 3C41: the fits for these galaxies indicate 
the presence of an unresolved nuclear component contribut- 
ing 37% and 24% of the total K-band flux respectively. It 
is interesting that these galaxies are the most luminous in 
our sample in the K-band, l ying furthest from the mean 
K— z relation. Rawlings et al. (1995) have already suggested 



that 3C22 is a reddened quasar observed close to the ra- 
dio galaxy / quasar divide, based upon the compactness of 
the infrared emission, the detection of broad H q em ission 
in their spectrum and that of Economou et al. (1995), and 



the very red colour of its spectral energy distribution. This 
hypothesis is supported by the fact that this source pos- 
sesses a bright radio core and a strong one-sided jet, fea- 
tures which are characteristic of quasars, but not generally 
of radio galaxies (see Paper I for further discussion). No such 
claim has yet been made for 3C41, although we note that 
it is one of the few sources in the sample to possess a clear 
radio jet (Paper I) . 
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Table 3. The percentage of point source component determined 
by the fitting routine for a given percentage of star added to 
the radio galaxy. This was carried out both for fitting where the 
characteristic radius was known from the HST data, and when it 
was a free parameter in the fit. 



Given the small size of these galaxies relative to the see- 
ing profile, to test the accuracy of the fitting procedure for 
determining the nuclear point source contribution, we added 
the emission from a star, at various percentages of the total 
flux density, to the infrared images of a number of the ra- 
dio galaxies. The resultant images were passed through our 
fitting procedure to determine the percentage nuclear point 
source recovered. The results are given in Table ^ and plot- 
ted in Figure ^ (crossed circles). For additions of a star of 
up to about 15% of the total flux density, the recovery of the 
point source component is almost total, giving confidence to 
the result that there is little or no requirement for a nuclear 
contribution in six of the galaxies in Figure ^ In particular, 
there is no evidence, either in the de Vaucouleurs fit or in 
the SED, for the strong (~ 45%) nuclear co ntrib ution to the 
emission of 3C65 suggested by Lacy et al. ( 1995 ). For larger 
levels of star added, typically only about 75% of the flux 



density is recovered as a point source. These results suggest 
that the true point source contribution for 3C41 will be ap- 
proximately 30%, and for 3C22 the best-fitting value will 
be about 50%, with an upper limit of as much as about 70% 
of the emission being associated with an unresolved nuclear 
component. The 'corrected' point source contributions are 
given in Table |^ together with the 90% confidence limits to 
this value as determined from the least-squares fits. 

Radial intensity profile fits can also be made for those 
galaxies which possess enhanced ultraviolet emission in the 
IfST images, since these too are dominated by the under- 
lying old stellar population at infrared wavelengths. 3C368 
and 3C247 were omitted from this analysis, the former be- 
cause of the presence of a foreground star within the galaxy 
envelope which cannot be accurately removed (see Paper I) , 
and the latter because it is surrounded by nearby compan- 
ions, too many to employ the quadrant method discussed 
above. The 5 highest redshift galaxies in the sample were 
also omitted due to low signal-to-noise ratios of the infrared 
emission at large radii. 

For the remaining 13 galaxies the infrared radial inten- 
sity profiles were measured, and the best-fitting matches 
were derived from the sum of an unresolved point source 
and a de Vaucouleurs profile. The characteristic radius of 
the de Vaucouleurs profile was allowed to vary to obtain the 
best match. The best-fitting models are compared with the 
observations in Figure ^, with the same notation as in Fig- 
ure ^, and the parameters of the models are given in Table ^ 
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Figure 6. The percentage of point source component determined 
by the fitting routine for a given percentage of star added to 
the radio galaxy. The crossed circles show the results when the 
characteristic radius was known from the HST data, and the filled 
triangles show the corresponding values for the galaxies for which 
the characteristic radius was a free parameter in the fit. 



Although the individual values of can only be determined 
to a typical accuracy of about 35%, the galaxy light profiles 
can again be well fitted using this model. 

The small size of the galaxies leads to some degeneracy 
between the characteristic radius and the point source con- 
tribution. We have repeated the test of adding a point source 
to the radio galaxies, allowing the characteristic radius to 
vary in the fit. The results are given in Table ^ and shown 
in Figure ^. The contribution of point source measured was 
lower than was added throughout, typically by about a fac- 
tor of two. What is happening is that the fitting routine is 
estimating a slightly smaller value for rc, and attributing the 
remainder of the point source flux to the galaxy profile. The 
error on the points in Figure ^ represents the scatter in the 
results for addition to different galaxies, and arises predom- 
inantly from the range of characteristic radii in the sample: 
the recovery of the point source is least effective for those 
sources with the smallest characteristic radii. However, even 
for small additions of a stellar component, some point source 
contribution was generally recovered, and therefore the very 
low point source percentages required for the galaxies in Fig- 
ure M (mean value < 4%) give confidence that any nuclear 
contribution must be small. The 'corrected' contributions 
of the unresolved emission to the K-band flux density are 
listed in Table |^ together with the 90% confidence limits 
on these values. Each of these galaxies is consistent with no 
point source component being present. 

For a significant fraction of the galaxies, the observed 
profiles become brighter than the predicted de Vaucouleurs 
profile at large radii, suggesting that at least some of these 
galaxies possess diffuse extended envelopes. The reality of 
these envelopes for any individual galaxy cannot be gauged 
because the signal-to-noise ratio of this extended emission 
is low and, more importantly, the errors at large radii are 
dominated by uncertainties in the background subtraction 
and so are not necessarily independent. We have therefore 
combined the radial intensity profiles of the galaxies. The 
radial profile of each galaxy was scaled to the same char- 
acteristic radius enabling the galaxy profiles to be summed 
directly. Only those galaxies for which 1.0 < < 2.0 were 
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Figure 7. Fits to tlie radial intensity profiles of the K— band images of the 13 SCR galaxies for which de Vaucouleurs profiles could 



not be determined from the 
notation is used as in Figure 



[ST observations. The characteristic radii have been determined from the best fitting profiles. The same 
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Figure 8. A combined radial intensity profile in the K-band for 
the 12 galaxies with 1.0 < ro < 2.0. Units on the y-axis are arbi- 
trary. The combined errors for each point are shown. The dotted 
line shows a combined de Vaucouleurs profile: halo emission is 
clearly visible at r > 2.5rc (ie. (r/rc)-'-''^ > 1.25). 



included: for galaxies with a smaller there remains a sig- 
nificant contribution due to the seeing at radii of 2 to 3 re, 
whilst for galaxies with re greater than 2, there is insufflcient 
signal at the largest radii to test for the presence of the halo. 
The scaled radial profiles of the remaining 12 galaxies were 
summed, weighting each galaxy equally, and the results are 
presented in Figure ^. 

At radii greater than about 2.4 re, that is, (r/rc)^^* ~ 
1.25, at which the effects of seeing should be small, a sig- 
nificant halo component has been detected. This envelope 
extends from a physical radius of typically about 35 kpc out 
to at least 60 kpc, and probably further. This extended en- 
velope may be dominated by the emission from only a small 
number of the radio galaxies, but it must be present in at 
least some of them. 

The halo emission is not an artefact associated with the 
fiat spectrum contribution to the K-band emission, for two 
reasons. Firstly, on the HST images (Paper 1), well over 95% 
of the fiat spectrum emission originates at radii less than 
35 kpc, the radius at which evidence for the halo component 
first appears; secondly, if the 12 galaxies are split into two 
groups according to the percentage of flat spectrum light in 
the K-band image, there is no significant difference between 
the strength of the halo component in the two samples. We 
have also studied the profiles of stars of the images to check 
that this halo is not an artefact associated with wings in the 
UKIRT point spread function. 

Such a halo component is typical of cD galaxies (e.g. 
Oemler 1976). Thuan and Rominishin (1981) compared the 
radial intensity profiles of brightest cluster galaxies in rich 
and poor clusters of galaxies, and showed that enhanced 
halo components were only seen around the central galaxies 
of rich clusters, and that they begin at radii between about 
25 and 40 kpc, comparable to the 35 kpc derived for the 3CR 
galaxies. Figure ^ strongly suggests that at least some of the 
3CR galaxies at redshifts z^l are cD galaxies in rich cluster 
environments. 



5 SURFACE PHOTOMETRY AND 

COMPARISON WITH LOW REDSHIFT 
ELLIPTICAL GALAXIES 

Elliptical galaxies at low redshift possess spectro- 
photometric properties which lie on the 'fundamental plane', 
in which their characteristic radii r^, average surface bright- 
nesses within their characteristic radius and velocity 
dispersions, ctv, are strongly correlated. The projection of 
the fundamental plane using only and {y)^, known as the 
Kormendy relation, also produces a tight correlation and 
avoids the need f or detailed spectroscopy ( Kormendy 1977 ; 
Schombert 1987; Oogorle & Hoessel 1991). This provides a 



powerful tool for investigating the stellar populations of high 
redshift elliptical galaxies. 

To compare the 3CR galaxies with existing data for 
elliptical galaxies at low redshift, the mean rest-frame V- 
band surface brightnesses within the de Vaucouleurs radii 
are required. Although in many cases an HST image at rest- 
frame V-band is available, these are dominated by the flat 
spectrum aligned component, and may also be influenced 
by dust extinction. As shown above, however, the K-band 
images are relatively unaffected except in a few cases; in the 
analysis that follows we omit 3C22 and 3C41 because of the 
presence of nuclear components which contribute to their 
K-band emission. 

We adopt the same assumptions as in Section ^, namely 
that the 3CR radio galaxies are passively evolving giant el- 
liptical galaxies which formed at a large redshift, with an ad- 
ditional fiat spectrum aligned component. We use the results 
of Section ^ to remove the contributions of any fiat spectrum 
or nuclear component from the K-band fiux densities of the 
galaxies. The K-band surface brightnesses thu s der ived are 
listed in Table |[ The Bruzual and Chariot ( |l993[ ) steUar 
synthesis codes were used to produce the SED of a passively 
evolving elliptical galaxy, which formed at Zi = 10, at the 
redshift of each 3CR source. These SEDs were then used to 
derive appropriate K-corrections and hence, taking account 
of the (1 + z)'* cosmological surface brightness dimming, the 
mean rest-frame V-band surface brightnesses. The varia- 
tion of the V-band flux density with age for each synthe- 
sised galaxy was determined, and was used to convert the 
derived V-band surface brightness into that which would be 
observed when the galaxy had continued to evolve passively 
through to zero redshift. This procedure resulted in surface 
brightnesses of the 3CR galaxies which can be compared di- 
rectly with low redshift giant ellipticals in the fundamental 
plane. 

The derived surface brightnesses are included in Ta- 
ble |2|, and are plotted against the de Vaucouleurs radius in 
Figure |^ (solid circles), together with data from low red- 
shift samples of brightest cluster galaxies in Abell clusters 
(crossed cir cles) and other ellipticals (crosses) taken from 
Schombert (1987). The radio galaxies lie along the funda- 
mental plane defined by the low redshift ellipticals, providing 
direct evidence that their stellar populations indeed formed 
at large redshifts and are passively evolving. The SEDs, and 
hence the K-corrections of either young or non-evolving stel- 
lar populations would be significantly different and the 3CR 
galaxies would not lie on the Kormendy relation. 

There is some evidence that the slope of the Kormendy 
relation defined by the 3CR galaxies is slightly steeper than 
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Figure 9. The V— band Kormendy relation for the SCR galaxies (solid circles) compared with low redshift ellipticals (crosses) and 
brightest cluster galaxies (crossed circles). The low redshift data are from Schombert (1987). The V-magnitude data of the SCR galaxies 
refer only to the stellar components, and have been derived from their K-magnitudes assuming that the stellar populations evolve 
passively. The ellipse indicates the error ellipse for the SCR galaxies for which characteristic radii were measured from the UKIRT data; 
those for which the measurements were obtained from the HST data have error ellipses less than half this size. The dashed line shows a 
line of constant total luminosity. 



that of the low redshift ellipticals. The dashed line on Fig- 
ure 1^ shows a line of constant total luminosity, that is, 
a line along which the product of the surface brightness 
and the square of the characteristic radius is constant. The 
SCR galaxies lie roughly along such a line suggesting that, 
if they all evolved passively until a redshift of zero, the re- 
sultant galaxies would have similar luminosities. This result 
is similar to that found by Lilly and Longair (1984) and the 
inclusion of the surface brightness information in Figure ^ 
confirms that the radio galaxies at z ~ 1 are old giant ellip- 
tical galaxies. It seems that the powerful SCR radio sources 
are only formed if their host galaxies attain a specific well- 
defined stellar mass. We return to this point in Section ^ 

The data in Figure ^ can be used to compare the range 
of values of of the high redshift SCR galaxies with those 
of the populations of low redshift brightest cluster galaxies 
and of other ellipticals. Normalised cumulative frequency 
distributions of for each of the three samples in Fig- 
ure]^ are plotted in Figure |lo| Two differences are apparent: 
firstly, the dispersion of Ve in the SCR galaxies is significantly 
smaller than that of the other two populations; secondly, the 
mean characteristic radius of the high redshift SCR galaxies 
{r^ = 14.6 ± 1.4) kpc is greater than that of the low red- 
shift ellipticals (rl — 8.2 ± 1.0) kpc, but about a factor of 
2.25 smaller than the low redshift brightest cluster galaxies 
(r%r = S2.7 ± 1.1) kpc. It should be noted that the character- 
istic sizes of the distant SCR radio galaxies are dependent 
upon the adopted value of go and, for go = 0, would be 
approximately 25% smaller. On the other hand, the small 
degeneracy between the detection of a point source compo- 
nent and the value of the characteristic radius, discussed in 
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Figure 10. Normalised cumulative frequency distributions of ro 
for the three different populations of galaxies: solid line - the high 
redshift SCR galaxies; dotted line - low redshift brightest cluster 
galaxies; dashed line - low redshift elliptical galaxies (excluding 
brightest cluster galaxies). The SCR galaxies have a much smaller 
spread in re than the low redshift elliptical galaxies. 

Section means that if there were to be a small nuclear 
contribution to these galaxies then the characteristic radii 
would be up to 20% larger than those derived. Neither of 
these effects would change the significance of the results. 

According to cannibalism models (e.g. Hausman and 
Ostriker 1978) the position of a galaxy along the fundamen- 
tal plane is interpreted as being related to its merger his- 
tory. Assuming that the mergers occur homologously (e.g. 
Schombert 1987 and references therein), when a small galaxy 
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Figure 11. A plot of the K-magnitudes measured for the 28 
SCR sources in the recent IRCAM3 observations against those of 
Lilly and Longair (1984). In the diagram, the new K magnitudes 
were reduced to the same apertures used by Lilly and Longair 
in their single element photometry. Most of the discrepancies are 
within the quoted errors on the original observations, but there 
are a few significant changes. 



is cannibalised by a much larger galaxy to produce a bound 
remnant, a large fraction of the kinetic energy of the merger 
is converted into binding energy of the remnant. The result 
of such a homologous merger is that the remnant has a larger 
radius and more diffuse morphology as compared with the 
original galaxy, thus moving the galaxy to th e right along 
the pr ojected fundamental plane in Figure ^ ( Hausman & 
Ostrike: 1978). In this respect, the larger size of the high 



redshift SCR galaxies as compared to low redshift elliptical 
galaxies is important: despite being observed at a much ear- 
lier cosmic epoch, the SCR galaxies at 2 ~ 1 appear to be 
highly evolved dynamically. We discuss the relative sizes of 
the distant SCR galaxies and the brightest cluster galaxies 
in Section ^ 

Within the SCR sample, it is interesting that SCS4 has 
the largest characteristic radius, suggesting that it has un- 
dergone most mergers, whilst SCS40 is the most compact 
object. The former of these is known to lie towards the cen- 



tre of a moderatel y rich cluster of galaxies (McCarthy 1988 



Best et al. 1997b), whilst there is pr actically no depolari 



sation of the radio lobes of the latter (Johnson et al. 1991:' 
suggesting that there is relatively little cluster gas surround- 
ing the galaxy. 



6 THE REVISED K-Z DIAGRAM 

Much evidence has been presented in the previous sections 
confirming the stellar nature of the infrared emission of the 
galaxies that host the SCR radio sources. A revised K— 2 
relation for the SCR galaxies has been derived, including a 
number of improveme nts ov er the original version presented 
by Lilly and Longair (1984): 



(i) The K-magnitudes measured by Lilly and Longair 
were obtained by single element aperture photometry using 
blind offsets and chopping. The new magnitudes show that, 
although in the majority of cases the 1984 infrared magni- 
tudes were remarkably accurate, in a small number of cases 



the on-source apertures were sufficiently large to contain 
emission from companions, or the off-source reference beam 
included a faint companion (Figure ^). The new imaging 
infrared observations allow these problems to be eliminated. 

(ii) Lilly and Longair had to make significant aperture 
corrections for many of the high redshift radio galaxies; 
in most cases, our magnitudes were measured through the 
adopted 9 arcsec diameter aperture. Any extended halos are 
generally only present at a radius r > S5kpc (4 arcsec at 
2 ~ 1), and so the contribution of any halo to the total 
K-band flux density within the 9 arcsec diameter aperture 
adopted for our photometry is unimportant. 

(iii) The sources SC22 and SC41 are not included with the 
narrow-line galaxies, but are classified with the broad line 
galaxies in which a nuclear component makes a significant 
contribution to the K-band fiux density. 

(iv) The contribution of the aligned component to the K- 
band magnitude has been removed by subtracting the ap- 
propriate fraction of flux density as estimated by our SED 
fltting in Section 3.1. The fraction of the total flux density 



that might be associated with any point source component 
has also been removed. The combination of these two cor- 
rections generally amounts to less than 20%, corresponding 
to about 0.2 magnitudes. The corrected Kcorr magnitudes 
are listed in Table |l|. 

For the radio galaxies with redsh ifts 2 < 0.6, the pho- 
tometry of Lilly and Longair ( 198^ ) was adopted. These 
galaxies are sufficiently bright that offsetting onto nearby 
faint companions would have had only a small effect upon 
the observed magnitudes, and there is no evidence for a 
strong alignment effect at low redshifts. The magnitudes 
of these galaxies have all been scaled to a 9 arcsec diam- 
eter aperture, using a radial light profile appropriate for low 
redshift radio galaxies (e.g. Lilly et al. 1984b): the correc- 
tions required are small, being less than 0.3 magnitudes in 
all cases, and generally less than 0.1 magnitudes. 

The revised K— 2 relation is shown in Figure ^ The 
broad-line radio galaxies, together with SC22 and SC41, 
are plotted as open circles. The predicted Vi—z tracks for 
various evolution models and cosmological parameters are 
also shown, normalised to match the low redshift (2 < 0.1) 
narrow-line radio galaxies. These tracks include a correction 
based upon a standard radial light curve for low redshift gi- 
ant elliptical galaxies, to account for the fact that the galax- 
ies are observed through an aperture of fixed angular size of 
9 arcsec, rather than fixed physical size. At redshifts 2 ~ 1, 
generally over 90% of the flux density is contained within 
this aperture. 

Figure 12k s hows a set of tracks for a go = universe, 
whilst Figure |12p shows the corresponding tracks in a qo = 
0.5 universe, in both cases the cosmological constant A being 
taken to be zero. In each case, the solid line represents the 
track of a non-evolving giant elliptical galaxy, obtained by 
redshifting the spectrum of a giant elliptical galaxy at the 
present epoch. The dashed, dotted, dash-dot lines show the 
predicted tracks of passively evolving galaxies formed in a 
1 Gyr burst at three different redshifts. A good fit to the 
data is obtained by the passively evolving models in the 
qo = 0.5 Universe. Little difference is found if the galaxies 
formed at redshift 2 = 5 or 2 = 20, but at redshifts 2 > 
1.2 the predicted track for galaxy formation at 2 = S is 
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Figure 12. Plots of K-magnitude vs redshift for the SCR radio galaxies. Broad line radio galaxies are marked by open circles as are 
3C22 and 3C41, and narrow line radio galaxies by filled squares. The error in the K— magnitude measurements is indicated by the vertical 
error bar associated with each point; in many cases, this is smaller than the plotted symbol. The upper plot corresponds to qo = and 
the lower plot to go = 0.5. The lines on each plot show the predicted K— z tracks of various evolutionary models and have been reduced 
to the magnitudes observed within a 9 arcsec diameter aperture at each redshift. The solid lines correspond to a non-evolving stellar 
population, whilst the dotted, dashed, and dash-dot lines are for passively evolving galaxies formed at Z[ =20, 5 and 3 respectively. The 
theoretical curves are normalised to match the low redshift (z < 0.1) narrow— line radio galaxies. 



significantly brighter tlian tlie observed data, indicating tiiat 
tiie stellar populations formed at an earlier cosmic epocii. 
None of tlie tracks in the go = plot provide a good fit to 
the observations. 



7 DISCUSSION 



7.1 The cluster environments of the distant SCR 
radio galaxies 

The amplitudes of the spatial cross-correlation functions 
for galaxies in the vicinity of near by SCR FRIT class ra- 
dio sources ( Fanaroff fc Riley 1974) are similar to t hose for 
normal elliptical galaxies ( [Prestage fc Peacock 1988 ) , whilst 
the host radio galaxies themselves have optical luminosi- 
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ties and characteristic sizes signific antly less than those o f 
first ranked Abell cluster galaxies (Lilly & Prestage f987). 



This indicates that low redshift FRIl sources typically lie 
in isolated environments or small groups. In contrast, mea- 



sures of the galaxy cross-correlation function ( Yates et al 



198S 
1991 

shifts. 



end an Abell clustering classification (Hill & Lilly 



) both indicate SCR galaxies with intermediate red- 
< 



0-3 z < 0.5, belong to environments which are 
about three times richer than those nearby. 

There is evidence that the high redshift SCR sources 
belong to moderate to rich cluster environments, in particu- 
lar: the detection of X-ray emission, in some cases resolved, 
from a number of these sources has been associated with a 



cooling flow in relatively dense intraclu s ter medium ( Craw- , 
ford & ^abian 1993|; [Worrall et al. 1994 [Crawford fc Fabian 



1995 |; iqrawford fc Fabian 199e| ; pickinson 1997D ; companion 
galaxies are seen around a number of the radio galaxie s in 
the narrow-line [Oil] S727 imaging by McCarthy ( 1988 ); in- 
frared imaging of distant radio galaxies provides large num- 
bers of cluster candidates, based upon selection by colour, 
and spectroscopic follow-up of individual SCR radio galaxies 
indeed shows them to lie in at least moderately rich clusters 
(e.g. Dickinson 1997); a surrounding medium of fairly high 
density is required to account for the large rotation measures 
of distant radio sources, and t o provide a working surface for 



Best et al 



their high radio luminosities (Carilli et al. 1997 
1997a|J 

In this paper we have provided further evidence in sup- 
port of this hypothesis. In Section ^ we showed that at least 
some of these galaxies possess halos characteristic of those 
seen around cD galaxies at low redshifts. Furthermore, the 
stellar masses of the SCR galaxies derived in Section pjl in- 
dicate that they are amongst the most massive galaxies at 
cosmic epochs at which they are observed. A comparison of 
the SCR galaxies with the brightest cluster galaxy (BCG) 
sample of Aragon-Salamanca et al. (199S) shows that the 



absolute K-magnitudes, and hence stellar masses, of the two 
samples at redshift z ^ 1 axe essentially the same. In con- 
trast, at low redshifts, BCGs are about a magnitude brighter 
than the SCR radio galaxies in absolute K-magnitude, con- 
sistent with the result that low redshift SCR galaxies lie in 
isolated environments of small groups. This difference is best 
illustrated by the fact that the K—z relation of the BCGs 



is consi stent with non- evolving stellar populations ( Aragon 
Salamaica et al. 199S). 

To understand the K—z relations, we must therefore 
unravel two 'cosmic conspiracies': 

(i) On the one hand, the tightness and slope of the SCR 
K—z relation suggests that the radio galaxies belong to a 
single population of objects which formed at large redshift 
and evolved passively over the redshift interval 2 > z > 0, 
such that a high redshift SCR galaxy would evolve into a low 
redshift SCR galaxy. This simple picture cannot be correct 
because we have argued that the galactic environments of 
FRII SCR radio galaxies change with redshift. High redshift 
SCR galaxies appear to lie in (proto-) cluster environments, 
whilst the nearby sources in the sample are found in isolated 
environments or small groups. 

(ii) On the other hand, the K—z relation of the BCGs 
suggests that their stellar populations are non-evolving, but 
this is clearly an unphysical picture. All passively evolving 
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Figure 13. A plot of the stellar mass derived by a fit to the 
broad band photometry of Section H against redshift for 26 SCR 



models for stellar populations suggest that they should be 
about a magnitude brighter in absolute K-magnitude at red- 
shift z = 1 as compared with z — 0. 

Since the stellar populations of the BCGs must have 
evolved in the same way as the SCR radio galaxies, the 
shape of their K—z relation reflects the fact that central clus- 
ter galaxies continue to accumulate matter through mergers 
and gas infall. Hierarchical clustering models for structure 
formation have suggested that the mass of BCGs can in- 
crease by a factor of up to 5 between a redshift of one and 
the present epoch (Aragon-Salamanca et al. 1997). In con- 
trast, the variation of the stellar masses of the SCR galaxies 
with redshift, as derived from the SED flts of Section S.l, 
shows no correlation (see Figure |l^) . 

In Section ^ we showed that the mean characteristic ra- 
dius of the distant SCR radio galaxies was approximately a 
factor o f 2.25 smaller than that of the nearby BCGs. Kor- 
mendy (1977) derived a relation between the integrated lu- 
minosity of a bright elliptical galaxy and its characteristic 
radius: Lint oc r^'^, or equivalently, Mstei oc {M/L)or°''^, 
where Mstei is the stellar mass of the galaxies and (M/L)o 
is the mass-to-light ratio of the stellar populations at red- 
shift zero. This result was confirmed for elliptical galaxies 
with Tc > 6kpc by Schombert (1987). This implies that the 
stellar masses of the SCR radio galaxies at redshift z ~ 1 
are only a factor of about 1.75 below those of nearby central 
galaxies in Abell clusters, assuming that their stellar popula- 
tions will have a similar mass-to-light ratio when they have 
evolved to redshift zero. When compared to the growth fac- 
tor of between S and 5 expected for BCGs over this redshift 
interval, it is clear that the distant SCR galaxies are mas- 
sive galaxies, even relative to the progenitors of the central 
galaxies of rich cluster. 



Intriguingly, West (1994) has shown that the comoving 



number density of powerful high redshift radio galaxies is of 
the same order of magnitude as that of rich Abell clusters, 
if account is taken of the fact that radio source lifetimes, 
~ 10* years, are much shorter than the age of the Universe. 
It is quite conceivable that all of the brightest galaxies in 
clusters today went through a phase of during which they 
were intense radio sources, and so must possess supermassive 
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black holes in their nuclei (e.g. van der Marel et al. 1997 and 
references therein). 

7.2 Understanding the SCR radio galaxies 

Since SCR radio galaxies at large redshifts lie at the centre 
of galaxy clusters, their masses would be expected to in- 
crease with cosmic epoch just like the BCGs; the 'passively 
evolving' K—z relation of these galaxies suggests, however, 
that this does not occur. This is a conspiracy: the galaxies 
sampled at high and low redshifts do not form a uniform 
population, as is indicated by the dramatic change in their 
galactic environments with redshift. To understand this ef- 
fect, we compare the K—z relation of the SCR galaxies with 
that of the only other complete sample of low-frequency 
selected radio sources, the 6C radio galaxies (Eales & Rawl- 
Eales etal. 199^ ). mportant differences 



ings 1996 



between the two samples, specifically: 

(i) At a given redshift, the SCR radio sources are about 
a factor of 5 more powerful than the 6C radio sources, since 
the latter are selected from a fainter radio survey. 

(ii) At low redshifts, z < 0.6, the SCR and 6C radio galax- 
ies have similar infrared luminosities. 

(iii) At redshifts 2^1, the (un-corrected) K-band lumi- 
nosities of the SCR galaxies are greater th an those of the 6C 
galaxies by, on average, 0.6 magnitudes (Eales et al. 1997) 



althougjli Lliyrt! is sigiiificaiiL overlap between Llie Lwo popu- 
lations. 



colours, imply much older ages for the stellar populations 
of the SCR radio galaxies than would be expected in this 
picture (eg. Stockton et al. 1995), and that it is difficult to 
assemble such massive, dynamically evolved systems as the 
SCR galaxies in a short cosmic time. 

It is likely, therefore, that both the 6C and SCR pop- 
ulations formed at large redshifts, and that the 6C galaxies 
are fainter than the SCR galaxies in the K-band at high 
redshift because they contain, on average, a lower mass of 
stars. The primary selection criteria for these samples of ra- 
dio galaxies are the limiting radio flux densities of the SCR 
and 6C catalogues, and so the different K—z relations of 
these two samples imply that the radio luminosities of the 
sources must depend upon the masses and environments of 
the host galaxies. 

According to conventional ideas about the origin of 
PR II radio sources, their luminosities are ultimately deter- 
mined by the powers of the beams of particles ejected by the 
AGN, and by the density of the surrounding interstellar and 
intergalactic gas which controls adiabatic radiation losses of 
the radio lobes. In turn, the intrinsic beam power of the 
source is determined by the mass of the central black hole, 
the mass of material accreting onto it, and the efficiency 
with which this is converted into beam energy. Given an 
abundant supply of fuel, it is a reasonable assumption that 
the total energy output of the central engine, whic h is dom- 
inated by the bulk kinetic power of the radio jets (Rawlings 



Saunders 1991), will be close to the Eddington limiting 



Points (ii) and (iii) reffect the fact that the K—z relation for 
6C radio galaxies is co nsistent with a non-evo lving popula- 
tion of radio galaxies ([Eales fc Rawlings 199^ ). 



Three possible reasons why the K-magnitudes of the 
SCR galaxies might be brighter than those of the 6C galaxies 
at z > 1 are: (i) there may be a strong direct or indirect AGN 
contrib ution to the K-band emission of the high redshift 



SCR gc laxies (e.g. Eales et al. 1997); (ii) the SCR galaxies 
may be younger than the 6C galaxies, and hence their stellar 
popula tions would be more luminous; (iii) the SCR galaxies 



luminosity. Rawlings and Saunders derived the total kinetic 
power of the jets in radio galaxies with redshifts z ^ 1 and 
found values corresponding to the Eddington limiting lumi- 
nosity of a black hole with a mass of a few times 10* Mq; 
Kaiser et al. (1997) suggest that the jet powers may be even 
higher. Theorists have long argued that the masses of the 
central black holes should be roughly proportional to the 
masse s of the host galaxies ( Soltan 1982|; [Efstathiou fc Rees 



198^; 3maU fc Blandford 1992; Haehnelt fc Rees 1993), and 



some evidence for such a correlation has recently been found 



for massive black holes in the nuclei of nearby galaxies (Ko- 



may be more massive than the 6C galaxies and thus contain rmendy fc Richstone 1991). Thus, the luminosity of a radio 



a greater number of stars. 

Concerning the first possibility, nearly 50% of the K- 
band emission of the SCR galaxies would have to be associ- 
ated with AGN activity. The results presented in this paper 
show that, with the exceptions of the two sources SC22 and 
SC41, the mean point source contribution to the total K- 
band flux density is only 4%, with secure upper limits of 
< 15%. We have discussed why line emission from gas ex- 



source is determined by three principal factors: the mass of 
the host galaxy, the availability of fuel, and the density of 
the environment. 

The gas which forms the fuel for radio sources may arise 
either through the infall of interstellar and intergalactic ma- 
terial towards the central regions of the galaxy, or as a re- 
s ult of an interaction o r merger with a gaseous dwarf galaxy 



(tHeckman et al. 1986t ^mith fc Heckman 198£| ; ^mith fc 



cited bj the AGN light is unlikely to be important, and have Heckman 199C| ; [Mcnon 1992| ). In high redshift proto-cluster 



estimated the mean contribution of the aligned emission to 
the K-band flux density to be about 10%. We conclude, 
therefore, that the sum of these effects will contribute at 
most about O.S of the 0.6 magnitudes difference between 
the two samples. 

In the second case, if the SCR galaxies were formed at 
a smaller redshift than the 6C galaxies, say Z{ ^ 2.5, their 
younger stellar populations would be more luminous than 
the corresponding 6C galaxies at redshifts z > 1, but by 
z ~ 0.6 they would have aged and the galaxies would no 
longer be significantly more luminous. The main objections 
to this picture are that features detected in the spectra of 
some of the galaxies, coupled with their very red infrared 



or young cluster environments, the central galaxies will be 
massive, and so contain amongst the most powerful central 
engines at that epoch. There will also be high merger rates 
and a plentiful supply of disturbed intracluster gas to fuel 
the central engin e and confine the radio lobes. Furthermore, 
EUingson et al. (1991) have found that the velocity disper- 



sions of galaxies around distant powerful radio sources are 
significantly lower (400 to 500 km s"'^) than those in com- 
parably rich clusters at low redshifts (500 to 1000 kms^^). 
These smaller relative velocities increase the efficiency of 
galaxy mergers. Thus, massive galaxies in young cluster en- 
vironments at high redshifts possess all the necessary ingre- 
dients for producing the most powerful radio sources. 
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Radio sources in galaxies of lower mass than the SCR 
galaxies are likely to be powered by less massive black holes 
and therefore to have lower intrinsic beam powers. It is rea- 
sonable to attribute the lower radio luminosity of the high 
redshift 6C radio sources as compared with the SCR sources 
to the fact that the 6C galax;ies are less massive and possess 
less massive black holes. The overlap between the two sam- 
ples on the K— 2; relation would then arise from the scatter in 
the correlation between the stellar mass of the galaxy and 
the mass of the black hole (e.g. Kormendy and Richstone 
1995); for example, a galaxy with a K-magnitude compa- 
rable to that of the SCR galaxies, but whose central black 
hole is less massive than predicted by the correlation, would 
host a 6C radio source. 

At redshifts z < 0.6, the lower radio luminosity of all 
radio sources means that the galaxies are no longer produc- 
ing radio beams with kinetic powers close to the Eddington 
limit. The analysis of Rawlings and Saunders (1991) con- 
firms that the beam powers of these sources are significantly 
lower than the values found in the radio galaxies at 2 ~ 1. 
Neither the masses of individual black holes, nor the galaxy 
clustering can decrease with cosmic epoch, and so the radio 
luminosities of low redshift sources must be limited not by 
the mass of the host galaxy, as are high redshift sources, 
but rather by the availability of gas in the host galaxy and 
its environment to fuel the central black hole. As a result, 
at low redshift a weaker correlation would be expected be- 
tween galaxy mass and radio luminosity. This may account 
for the similarity of the K-band magnitudes of the SCR and 
6C galaxies at these redshifts. 

This decrease in the fuelling gas density with cosmic 
epoch has been argued on different grounds by other authors 
(e.g. Rees 1990 and references therein). It is also consistent 
with recent results on the evolution of the fraction of baryons 
in the form of neutral hydrogen, and the enrichment of the 
heavy elements, which indicate that the quantity of inter- 
stellar and intracluster materi al has decreased si gnificantly 
from redshifts 2; ~ 2 to z = (Madau et al. 1996). 



tion rates a radiatively efficient thin accretion disk forms. 
Studying the core spectra of the radio lo ud qu asar SC27S 
and the FRI source M87, Reynolds et al. (1996) argue that 



7.3 



Cosmic evolution of the radio source 
population 



We have argued that radio galaxies at high redshifts lie in 
young cluster environments, and that such environments 
possess the ideal ingredients for producing powerful radio 
sources. Why then are the most powerful FR II radio sources 
at the present epoch not also associated with central clus- 
ter galaxies? The radio sources hosted by nearby BCGs are 
almost invariably 'edge-darkened' FRI type radio sources. 
FRIs do not possess large luminous radio lobes; instead, 
the jets are decelerated significantly in the inner kiloparsec, 
probably due to entrainment of external material (Laing 
199S and references therein), implying a relatively low ki- 
netic power for FRI jets. How do the environments of low 
redshift central cluster galaxies restrict, in some way, the 
central engines of these galaxies to producing only the lower 
power FR I jets? 



Narayan and Yi (1995) have shown that if the accretion 



rate on to a black hole falls below a critical value given by 
Merit ~ Q^MEdd (where a is the viscosity parameter which 
has a value of order 0.1) then accretion on to the black hole 
occurs in an advection dominated mode. At higher accre- 



the difference between FRI and FRII sources may be that 
the central engines of FR I sources accrete matter at a rate 
below the critical value. The large difference in jet luminos- 
ity between FRI and FRII type sources would then arise 
from a much smaller difference in the mass accretion rate, 
coupled with a ch ange i n the accretion mechanism. 

Baum et al. ( 1992| ) showed that the gas surrounding 
FRI radio galaxies is of low angular momentum, with no 
evidence for recent mergers, whilst that surrounding FRII 
sources is generally of higher angular momentum and fre- 
quently exhibits tidal tails or other evidence for mergers. 
Rich regular clusters at low redshifts have had time to viri- 
alise the spatial distribution of the galaxies, and for the 
intracluster gas to take up an equilibrium configuration 
within the cluster grav itational potential. Studi es of the 
Butcher-Oemler effect ( [Butcher fc Oemler 197^) indicate 



that, compared to high redshift clusters, there is a dearth of 
gas-rich galaxies close to the centre of low redshift clusters 
which might merge with, and fuel, the central galaxy. Fur- 
thermore, the high velocity dispersion of the galaxies greatly 
reduces the merger efficiency. The evolution and virialisa- 
tion of low redshift clusters is therefore likely to result in 
a more restricted fuelling supply for the central engines of 
BCGs, and the fuelling gas may arise solely through steady 
accretion from the interstellar or intracluster media (see also 
Baum et al. 1995). The most powerful FRII radio sources at 
small redshifts are generally found instead in small groups 
of galaxies: in these groups, velocity dispersions are small, 
mergers relatively frequent, and gas-rich galaxies plentiful. 
Occasionally, dramatic mergers will occur in low redshift 
BCGs, producing radio sources of comparable radio power 
to the z 1 SCR sources, and we ascribe exceptional sources 
such as Cygnus A to this process (see also Barthel and Ar- 
naud, 1996). 

These considerations suggest an explanation for the 
constancy of the stellar masses of the SCR galaxies through- 
out the redshift range O.OS < z < 1.8. This mass may be in- 
terpreted as the typical mass which a giant elliptical galaxy 
attains before its galactic and gaseous environment is viri- 
alised. Galaxies of mass greater than a few times 10^^ 
will have undergone greater dynamical evolution, and gener- 
ally live in virialised environments where the reduced supply 
of fuelling gas to the central regions results in FRI sources 
being formed. This leads naturally to a correlation between 
the redshift of a SCR FRII galaxy and the richness of its 
environment: galaxies in highly overdense environments ac- 
cumulate matter the fastest, and therefore reach this 'FRII 
upper mass limit' at early cosmic epochs, whilst those which 
lie in less dense environments evolve more slowly and can 
form powerful FRII radio galaxies at the current epoch. 

In the model presented here, the location of the dis- 
tant SCR galaxies in the low redshift fundamental plane 
can be understood. All galaxies evolve through accumula- 
tion of matter and move to the right along the fundamental 
plane. For a galaxy to form a radio source at 2: ~ 1 of suffi- 
cient radio power to be in the SCR sample, the galaxy (and 
hence its central black hole) must have attained a certain 
mass, and therefore a certain luminosity, which is close to 
the dashed line on Figure The small scatter in their char- 
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acteristic radii follows naturally from the galaxies being seen 
at a similar point in their evolutionary history. 

This model is also fully consistent with the redshift evo- 
lution of the alignment effect seen in powerful radio galax- 
ies (e.g. McCarthy 1993, Paper III). At high redshifts, the 
sources lie in environments in which the surrounding gas is 
unsettled, distributed throughout the forming cluster, and of 
relatively high density. Most models of the alignment effect 
depend critically upon the availability of cool gas and so, 
at intermediate and low redshifts where the environments 
of the sources have much lower gas densities and the gas is 
more settled, there is a decrease in the luminosity of aligned 
emission. 

These results suggest that the highest redshift (z > 2) 
radio galaxies, like the SCR galaxies, are amongst the most 
massive systems at their observed epoch (see also Pentericci 
et al. 1997). These galaxies must have accumulated their 
mass at early cosmic epoch, and so must belong to envi- 
ronments in which the most massive systems are formed by 
a redshift of zero. The dense surrounding environment and 
the abundant supply of fuel required are consistent with 
the detection of large Lya halos around these objects (e.g. 
Rottgering and Miley 1996 and references therein), and the 
very large de polarisation and ro tation measures of their ra- 
dio emission (Carilli et al. 1997). 

A combination of luminosity evolution of the radio 
source population, and the greater availability of gas at 
high redshifts increasing the probability of a powerful ra- 
dio source being formed, provides a natural explanation for 
the increase in the comoving number density of radio galax- 
ies and quasars out to redshifts z ~ 2. The levelling off 
and decline in the comoving number density of powerful ra- 
dio sources beyond this redshifts z > 2 can plausibly be 
explained in terms of the fact that systems which are rich 
enough to produce sufficiently massive host galaxies, and 
consequently massive enough black holes in their nuclei, had 
not had time to form by that epoch. 

At the same time, these results also provide an expla- 
nation of why other populations of active galaxies, such as 
the optically selected quasars and the bright X-ray galaxies, 
show the same evoluti onary behaviour as the powerful radio 
galaxies ( Dunlop 1994 ) : since the maximum radio luminosity 
is governed by Eddington limited accretion, this same pro- 
cess will define the limiting optical and X-ray luminosities. 
Therefore, the cosmic evolution of all classes of active galaxy 
are governed by the same two parameters: the evolution of 
black holes masses, and the decrease in the availability of 
fuelling gas with cosmic epoch. 



• Except in two cases (3C22 and 3C41), there is no evi- 
dence for a nuclear point source contribution to the K-band 
emission at > 15% of the total flux density at that wave- 
length. 

• At least some of the high redshift 3CR galaxies possess 
extended envelopes, similar to the halos seen around nearby 
cD galaxies. Such envelopes are not seen around low redshift 
3CR galaxies, indicating a difference between the high and 
low redshift populations. The high redshift 3CR galaxies 
appear to live in environments which will form rich clusters 
by a redshift of zero. 

• The 3CR radio galaxies have a similar stellar mass over 
the complete redshift range of the sample. The less powerful 
6C radio sources have lower stellar masses at high redshifts; 
the black holes associated with these galaxies are likely to 
be correspondingly less massive, producing proportionally 
weaker radio sources. 

• Powerful radio sources are not observed in the most 
massive galaxies at low redshift. A reduced supply of fu- 
elling gas for the central engines, possibly associated with 
the lower merger rates and lack of gas-rich galaxies towards 
the centre of virialised clusters, results in lower luminosity 
FRI type radio sources being preferentially formed. 

• If the 3CR galaxies evolve in the same way as brightest 
cluster galaxies, then the evolution in the comoving number 
density of powerful radio sources can be naturally under- 
stood: out to redshifts 2; ~ 2, the comoving number density 
of powerful radio sources increases through a combination of 
luminosity and density evolution associated with the greater 
mass accretion rates and the greater probability of a strong 
radio source being formed; at redshifts z > 2, the levelling 
off and decline of the comoving number density would be due 
to negative density evolution, associated with the decreas- 
ing number of galaxies which lie in environments sufficiently 
rich enough for them to have accumulated enough matter to 
produce a powerful radio source by that epoch. 

• The commonly-used 'uniform population, closed box' 
galaxy evolution models are not appropriate for interpre- 
tations of the K— 2; relations. Powerful radio galaxies se- 
lected at high and low redshifts have different evolution- 
ary histories, but must contain a similar mass of stars, a 
few times lO^'^M©, and so conspire to produce the 'pas- 
sively evolving' K— z relation observed. The apparently sim- 
ple 'no-evolution' shape of the K— z relation for brightest 
cluster galaxies and lower power radio galaxies is the result 
of a further conspiracy in which these galaxies become in- 
trinsically fainter due to stellar evolution but also brighten 
because of mass increase through mergers. The shapes of 
these relations can be used to provide information about 
the merger histories of massive galaxies in clusters. 



8 CONCLUSIONS 

The main conclusions of the study of the old stellar popula- 
tions of these 3CR galaxies are as follows: 

• The magnitudes, colours and location of the distant 
3CR galaxies on the projected fundamental plane for ellip- 
tical galaxies indicate that their stellar populations formed 
at large redshift and are passively evolving. Their spectral 
energy distributions can be well matched using such an old 
stellar population together with an aligned component of 
relatively flat spectrum. 
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